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ABSTRACT 

In order to examine the relative importance of powerful starbursts and 
Compton-thick AGNs in NGC 6240, we have obtained mid-infrared images and 
low-resolution spectra of the galaxy with sub-arcsecond spatial resolution using 
the Keck Telescopes. Despite the high spatial resolution (~ 200 pc) of our data, 
no signature of the hidden AGNs has been detected in the mid-infrared. The 
southern nucleus, which we show provides 80-90% of the total 8-25 /xm luminos- 
ity of the system, has a mid-infrared spectrum and a mid-/far-infrared spectral 
energy distribution consistent with starbursts. At the same time, however, it is 
also possible to attribute up to 60% of the bolometric luminosity to an AGN, con- 
sistent with X-ray observations, if the AGN is heavily obscured and emits mostly 
in the far-infrared. This ambiguity arises because the intrinsic variation of prop- 
erties among a given galaxy population (e.g., starbursts) introduces at least a 
factor of a few uncertainty even into the most robust AGN-starburst diagnos- 
tics. We conclude that with present observations it is not possible to determine 
the dominant power source in galaxies when AGN and starburst luminosities are 
within a factor of a few of each other. 
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1. INTRODUCTION 

NGC 6240 is a nearby (D ~ 100 Mpc)^ system of interacting galaxies which epitomizes 
the AGN-starburst controversy surrounding ultra-luminous infrared galaxies (ULIRGs)^. To- 
gether with Arp 220 (Soifer et al. 1984), NGC 6240 was one of the first few galaxies discovered 
by IRAS to have an extremely large (> 10^^-10^'^ Lq) and dominant (> 0.9 Lhoi) infrared lu- 
minosity (Wright, Joseph, & Meikle 1984). Its highly disturbed morphology (Fosbury & Wall 
1979) and two nuclei separated by less than 2" (< 1 kpc) (Pried & Schulz 1983) indicate that 
this is a system of two merging galaxies. Based on the shape of the infrared spectral energy 
distribution (SED) and the large infrared size (> 3 kpc) inferred initially, it was argued that 
powerful starbursts induced by galaxy interaction generate the large infrared luminosity by 
heating interstellar dust with the UV radiation from young massive stars (Wright, Joseph, 
& Meikle 1984; Joseph & Wright 1985). 

Although the size argument later turned out to be incorrect when the mid-infrared size 
was measured to be < 500 pc by Wynn- Williams & Becklin (1993), this starburst hypothesis 
for NGC 6240 has subsequently gained wide support from a series of follow-up observations. 
The most crucial pieces of evidence are, (1) the deep CO absorption bandheads, (2) the 
strong PAH features, and (3) the starburst-like radio-infrared luminosity ratio. The deep 
CO absorption bandheads detected in the near- infrared clearly indicate that the near-infrared 
continuum of NGC 6240 is dominated by starlight (Rieke et al. 1985; Lester, Harvey, & Carr 
1988; Ridgway, Wynn- Williams, & Beckhn 1994; Shier, Rieke, & Rieke 1996; Sugai et al. 
1997; Tccza et al. 2000). The most recent study by Tecza et al. (2000) shows that these 
CO absorption bandheads are likely to be produced by late K or early M supergiants, and 
that the starburst population associated with these supergiants, estimated to be triggered 
~ 2 X 10^ years ago and lasting ~ 5 x 10^ years, are powerful enough to provide 30-100% 
of the bolometric luminosity. The strong PAH features seen in the mid-infrared spectra of 
NGC 6240 (Smith, Aitken, & Roche 1989; Genzel et al. 1998; Dudley 1999; Rigopoulou et 
al. 1999; Imanishi & Dudley 2000; Lutz et al. 2003) also suggest starbursts as the main 
luminosity source since they are much weaker in Active Galactic Nuclei (AGN) galaxies. 



^We adopt a distance of 103.86 Mpc derived by Sanders et al. (2003) based on a redshift of 0.0243 
(Solomon et al. 1997) and the cosmic attractor model of Mould et al. (2000). The assumed cosmological 
parameters are Qm = 0.3, = 0.7, and Hq = 75 km Mpc~^, respectively. At this distance and redshift, 
1" subtends 490 pc. 

^With an infrared luminosity Ljr (= L(8-1000 /zm) as defined by Sanders & Mirabel (1996)) of 7.1 x 
1O"L0 (Sanders et al. 2003), NGC 6240 is a luminous infrared galaxy (LIRG; Lm > lO^^L©) rather than 
an ultra-luminous infrared galaxy (ULIRG; Lm > in the strict sense. However, we do not make 

this distinction here because NGC 6240 shares many of the ULIRG properties. 
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Finally, the radio-infrared luminosity ratio of NGC 6240, averaged over the whole galaxy, is 
also starburst-like (Colbert, Wilson, & Bland-Hawthorn 1994; Tecza et al. 2000) although 
the existence of compact radio cores (< 26 pc) in the two nuclei may suggest some AGN 
radio emission (Colbert, Wilson, & Bland-Hawthorn 1994; Beswick et al. 2001; Gallimore 
& Beswick 2004). Together with the starburst-driven superwind seen both in the visual 
(Heckman, Armus, & Miley 1987; Armus, Heckman, & Miley 1990; Heckman, Armus, & 
Miley 1990) and in the soft X-ray (Schulz et al. 1998; Komossa, Schulz, & Greiner 1998; 
Iwasawa & Comastri 1998; Lira et al. 2002), the evidence for the powerful starbursts is 
indeed strong. 

This starburst picture was seriously challenged by a series of the strong hard X-ray 
detections, which suggest the existence of a QSO-like AGN in NGC 6240. An initial hint 
came from the ASCA hard X-ray spectrum, which showed an extremely flat continuum above 
3 keV with a strong iron K hne complex (Mitsuda 1995; Kii et al. 1997; Iwasawa & Comastri 
1998; Nakagawa et al. 1999). This suggests that what we are seeing is reflected X-ray emission 
from an AGN obscured by a Compton-thick [Nh > 2 x lO^^cm^^) material, and that the 
intrinsic luminosity of this AGN could be as high as 2.6 x 10^^ L0(1O'^^ erg s^^) (Iwasawa & 
Comastri 1998). This conclusion was later confirmed by the Beppo-SAX observations, which 
showed that the X-ray spectrum of NGC 6240 clearly exhibits the emergence of a powerful 
hard X-ray continuum above 10 keV from the obscured AGN (Vignati et al. 1999). If we 
take the estimated intrinsic 2-10 keV nuclear luminosity of > 10^^ erg s~^ (absorption- 
corrected based on the Beppo-SAX spectrum extending up to 100 keV) and assume that it 
is less than 10% of the AGN bolometric luminosity (-^2-ioitey/-^6oZ ~ 0.03 — 0.1 for Seyferts 
and QSOs according to Iwasawa et al. (2001)), this sets a lower limit of 10^^ erg s~^ on 
the AGN luminosity in NGC 6240, which is consistent with the ASCA result. Given various 
uncertainties, the AGN can also produce 50-100% of the bolometric luminosity of NGC 6240. 
Recently, Chandra hard X-ray images of NGC 6240 have shown that this hard X-ray emission 
originates from two nuclei (Komossa et al. 2003). A detailed X-ray spectroscopic study has 
also been performed by Boiler et al. (2003) and Netzer et al. (2005) based on the XMM- 
Newton X-ray data. 

In the last few years, NGC 6240 has been studied with increasingly high spatial resolu- 
tion. Gerssen et al. (2004) resolved both nuclei into separate components using HST obser- 
vations at visual wavelengths, but they find no clear sign of the two AGNs. They emphasize, 
however, the importance of X-ray data for identifying AGNs in highly dust-enshrouded envi- 
ronments. Very Long Basehne Array (VLBA) observations at radio frequencies reveal three 
compact sources in NGC 6240 (Gallimore & Beswick 2004), two of which are associated with 
the X-ray nuclei. 
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The two conflicting pictures put NGC 6240 at the heart of the AGN-starburst contro- 
versy within ULIRGs: that is, do AGNs or starbursts dominate their luminosity output? In 
recent years, the prevaihng view has been that it is starbursts that dominate. The majority 
of ULIRGs shows no obvious sign of a hidden AGN, and the observed starburst components 
seem powerful enough to provide the bolometric luminosity of these galaxies when corrected 
for extinction (Genzel et al. 1998). However, NGC 6240 is a notable anomaly in this pic- 
ture with strong evidence for a dust-obscured powerful AGN. If this AGN indeed provides a 
substantial fraction of the infrared luminosity, the validity of various AGN-starbursts power 
diagnostics must be re-examined since many of these diagnostics indicate NGC 6240 to be 
a typical starburst galaxy. Most recently, the existence of a Compton-thick AGN was also 
hinted for Arp 220, another ULIRG with no sign of AGN signature in the infrared, through 
the detection of the Fc K emission in the hard X-ray (Iwasawa et al. 2005). The SED mod- 
eling of Arp 220 by Spoon ct al. (2004) suggests that such a heavily obscured AGN might 
contribute significantly to the bolometric luminosity. 

Here we present mid-infrared images and spectra of NGC 6240 obtained with the Keck 
Telescopes with sub-arcsecond spatial resolutions. The Keck Telescopes routinely deliver a 
diffraction-limited spatial resolution of 0'.'3-0'.'5 in the mid-infrared, which corresponds to 
150-250 pc at the distance of NGC 6240. Such a high spatial resolution is essential to probe 
a complicated merging system like NGC 6240, in which the two nuclei are separated by only 
l'/4-l'.'5 (l'/5 in VLB A — (Gallimore & Beswick 2004); 1'.'4 in X-rays — (Komossa et al. 
2003)). 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Keck II/MIRLIN Mid-Infrared Images 

The mid-infrared images of NGC 6240 were taken on UT 1998 March 19 with the 
MIRLIN camera (Ressler et al. 1994) on the Keck II 10m Telescope on Mauna Kca in 
Hawaii. The camera uses a 128x128 Si:As BIB array, and was attached to the f/40 bent 
Cassegrain visitor port, producing a pixel scale of O'.'ISS pixel"^ with a field of view 18" 
on a side. Secondary-mirror chopping and telescope nodding were employed for effective 
subtraction of the sky background and instrumental noise. The observational methods and 
characteristics as well as the photometric calibration and data reduction are the same as 
those used by Soifer et al. (1999, 2000) and are described therein. 

Observations of NGC 6240 were made in the seven mid-infrared filters listed in Soifer 
et al. (1999), where the central wavelengths and filter full widths are given. The central 
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wavelengths are also listed in Table 1 together with the on-source integration times. The 
mean airmass was 1.07. The MIRLIN observations emphasized photometry rather than spa- 
tial structure, and as a result no strict size limit can be placed on the observations of the 
nucleus of NGC 6240. In particular, since no star was observed close in time, we do not 
have accurate estimates of the point spread function (PSF) and seeing. However, the bright 
compact southern nucleus of NGC 6240 (Figure 1) has FWHMs of 0'.'5-0'.'7 at 7.9-24.5 /xm, 
indicating that the intrinsic seeing was at least this good. 

2.2. Keck I/LWS Mid-Infrared Spectra 

The low-resolution (AA/A ~ 50) mid-infrared spectra of NGC 6240 were obtained with 
the Long Wavelength Spectrometer (LWS; Jones & Puetter (1993)) on the Keck 1 10m 
Telescope. LWS uses a 128x128 BIB array, and is attached to the f/25 forward Cassegrain 
focus of the telescope, producing a pixel scale of O'.'OS pixel"^ with a field of view 10" on a side. 
The sht width was 0^48 (6 pixels). The observing methods, instrumental characteristics, and 
data reduction used to obtain the mid-infrared spectra are the same as were used by Soifer 
et al. (2002), and are described therein. 

The spectra of NGC 6240 covering 8.8 //m-13.1 pim were taken on UT 2000 May 21 while 
those covering 7.5 /xm-12.1 /xm were taken on UT 2000 May 22. The total integration times 
are 147 minutes and 59 minutes, respectively. The sht position angle was 16° (east of north), 
which is along the direction connecting the two nuclei. The NGC 6240 spectra were divided 
by the spectra of HR 5340 and HR 6406, and multiplied by a blackbody spectrum at the 
star's effective temperature. For wavelength calibration, we adopted the spectral dispersion 
of 0.0375 fim pixel^^ derived from our other LWS data, and adjusted the offset by using the 
PAH 11.3 //.m feature seen redshifted to 11.6 /^.m in the NGC 6240 spectra. 

3. RESULTS 

3.1. Mid-Infrcired Morphology and Photometry 

Figure 1 shows the mid-infrared images. The images clearly show two compact nuclei 
separated by ~ l'/6 ± O'/l along the position angle of 16° (east of north). The southern 
nucleus dominates the mid- infrared luminosity output. It appears round and symmetric 
although some north-south elongation is seen in the 11.7 and 12.5 /im images, which are of 
the highest signal-to-noise ratio. The sub-components of the two main sources seen by Max 
et al. (2005) in the near infrared and by Gallimore & Beswick (2004) with the VLBA cannot 
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be accurately resolved with the present mid-infrared resolution, but the elongations in the 
mid-infrared images hint at this structure. 

Ignoring the elongation, we have estimated the source size of the southern nucleus by 
fitting a circularly symmetric Gaussian to the observations. Since we do not have the data 
to determine the seeing at the time of the NGC 6240 observations, we simply subtracted, 
in quadrature, the diffraction-limited beam sizes from the measured Gaussian FWHMs and 
derived the most indicative upper limits on the intrinsic size. There was, however, no evidence 
that its size at any of the mid-infrared wavelengths exceeded 200 pc (0'.'4) after diffraction 
broadening of the images is accounted for. 

Table 1 lists the continuum flux densities of the two nuclei measured in a circular 1" 
diameter beam as well as the total flux densities for the whole system measured in a 4" 
diameter circular beam. The total flux densities at 12.5 and 24.5 /im measured by MIRLIN 
are 0.5 and 3.3 Jy, respectively, which are close to the IRAS flux densities of 0.6 and 3.4 Jy 
at 12 and 25 //m. This indicates that these two compact nuclei provide almost all the mid- 
infrared luminosity of NGC 6240. The southern nucleus is especially luminous, providing ~ 
80-90% of the total flux at these wavelengths. 

Table 2 compares the flux ratio of the two nuclei at a number of wavelengths from 
X-ray to radio. The flux ratio increases by a factor of two from near-infrared to mid- 
infrared, indicating that the southern nucleus is slightly redderr than the northern one in 
this wavelength range. However, the ratios are quite similar in the X-ray and radio (~2-3). 

Figure 2 compares the MIRLIN 12.5 /xm image, which has the highest signal-to-noise 
ratio, with the HST/NICMOS 1.1 and 2.2 //m images (Scoville et al. 2000). The positions 
of the two radio nuclei (GaUimore & Beswick 2004) are also overlaid (crosses) with the 
assumption that the positions of the southern nucleus are spatially coincident. Note that it 
is not possible to perform independent astrometric calibration with the mid-infrared images 
because NGC 6240 is the only detected source in the small held of view of the MIRLIN 
images. As the flgure shows, the separation and position angle of the two nuclei at 12.5 /im 
are almost exactly the same as those at radio, and hence X-ray, wavelengths. This suggests 
that we are looking at the same components in the mid-infrared as we see at the other 
wavelengths. 

The HST/NICMOS images, on the other hand, show a slightly larger separation be- 
tween the nuclei, and the northern nucleus is signiflcantly displaced from the radio position. 
In fact, Gerssen et al. (2004) have identifled a secondary peak in the HST/NICMOS im- 
ages of the northern nucleus of NGC 6240 close to the position of the radio source. Its 
faintness in the 1.1 iim image (Figure 2f) indicates that this source is likely to be heav- 
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ily extincted by dust absorption in the galaxy. The positional coincidence argues that this 
fainter secondary component is responsible for the mid-infrared and radio luminosity of the 
northern nucleus. The southern nucleus also displays a substantial north-south elongation in 
the HST/NICMOS images due to another secondary component seen conspicuously in the 
1.1 iim image (Figure 2f). 

3.2. Mid-Infrared Spectra 

Figure 3a shows the LWS two-dimensional spectrum (wavelength vs. location along the 
slit) of NGC 6240 covering the wavelength range from 9.9 /xm to 12.4 /xm in the restframe. 
The two spectra have been combined to form this two-dimensional spectrum. The actual 
slit position is shown in Figure 2. The southern nucleus is clearly visible, but the northern 
nucleus at the position of the upper white dashed line is too faint to see. We therefore 
exclude the northern nucleus from the following discussion. The southern nucleus has a 
detectable but faint continuum at < 11 /im, which rises sharply toward longer wavelengths. 
The prominent bright spot between 11 and 11.5 iim is due to the 11.3 iim PAH feature. 

Figure 3b shows the same spectral image with the continuum subtracted from the south- 
ern nucleus. The continuum subtraction was done in the manner described in Soifer et al. 
(2002). The 11.3 /xm PAH feature is clearly seen. However, the line emitting region seems 
to be displaced from the continuum by 0'/2 (~90 pc) downward, which corresponds to the 
south-west direction. This can be clearly seen when the spatial profiles of the continuum 
and PAH feature are compared (Figure 3c). Incidentally, this separation is similar to that 
of the two near-infrared peaks in the southern nucleus seen in the HST/NICMOS 1.1 
image (Figure 2e and f), suggesting a possible connection between the PAH peak and the 
secondary near-infrared peak. We also note that there is no obvious detection of the PAH 
emission from the molecular gas concentration between the two nuclei detected with the 
near-infrared H2 emission (Herbst et al. 1990; van der Werf et al. 1993; Sugai et al. 1997; 
Tecza et al. 2000; Bogdanovic et al. 2003) and with the millimeter CO emission (Tacconi et 
al. 1999). 

Figure 4a shows the LWS spectrum of the southern nucleus over the full wavelength 
range from 7.3 to 12.8 /im in the rest frame, together with the photometric measurements 
with MIRLIN. To produce the one dimensional spectrum, the signals within 7 pixels (0'.'56) 
from the continuum center were coadded. This extraction aperture is large enough to contain 
most of the continuum and PAH emission, and therefore the 0'.'2 offset of the PAH emitting 
region mentioned above docs not affect the resultant spectrum significantly (the effect on 
the PAH flux is a ~ 5% reduction at most). The flux level of the spectrum was adjusted 
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such that it makes a reasonable match to the broad band MIRLIN photometric points. A 
strong PAH 11.3 /xm feature is seen, and its hne flux is measured to be 8.6 x 10~^^ W m~^. 
Also, the sharp drop of flux at the short-wavelength end of the spectrum suggests that we 
are seeing the peak of the 7.7 fim PAH feature while the shoulder to the right indicates the 
existence of the 8.6 /xm PAH feature. 

For comparison, the ISOCAM CVF spectrum presented by Lutz et al. (2003) is shown in 
Figure 4b. (Lutz et al. (2003) also give ISOSWS data, but these are high resolution spectra 
targeted on individual lines). Since the ISOCAM CVF spectrum contains the fluxes from 
both nuclei, it is compared with the total flux measurements by MIRLIN with a 4" diameter 
circular beam. The two data sets show good agreement, confirming the accuracy of the 
photometric cahbration. The LWS spectrum of the southern nucleus is over-plotted, showing 
that the overall continuum shape (e.g., the depth of the 10 /im trough) is quite similar to 
that of the ISOCAM spectrum although the emission features (e.g., 11.3 //m PAH) are much 
weaker. 



4. DISCUSSION 

4.1. The Southern Nucleus 

As is the case with other ULIRGs (Soifer et al. 2000), the southern nucleus of NGC 6240, 
which dominates the luminosity output, is characterized by an extremely high surface bright- 
ness. If we take the intrinsic source area to be a circular disk with a diameter of 200 pc 
on the sky, the bolometric luminosity of 7 x 10^^ Lq translates into a surface brightness 
of 2 X 10^^ Lq kpc^^. This is at the high end of surface brightnesses of typical starburst 
galaxies (Soifer et al. 2000; Evans et al. 2003). The real surface brightness is probably even 
higher since we have taken an upper limit on the size of the continuum emitting region. 

Despite this high surface brightness, the mid-infrared spectrum of the southern nucleus 
is similar to those of the much less luminous dusty starburst galaxies such as M 82 and 
NGC 253. Figure 5 compares the LWS spectrum of the southern nucleus with the ISOSWS 
spectra of M 82 and NGC 253 (Sturm et al. 2000). These spectra can be compared directly 
since they both sample the light within an area of a few hundred parsec in diameter around 
the nucleus. The match is reasonable in terms of the 10 fira trough depth and PAH feature 
strengths. On the basis of the similarity of the mid-infrared spectrum to those of M 82 and 
NGC 253, we would argue that the mid-infrared emission originates from starbursts. 

The good agreement between the LWS and ISOCAM spectra seen in Figure 5 can be 
interpreted as a consequence of the southern nucleus dominating the total luminosity of 
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the system. For this reason, we will not distinguish between the LWS spectrum and the 
ISOCAM spectrum in the following discussion once an appropriate scaling is applied to the 
latter to match the former. 

Although the 7-13 //m spectra match well, the spectrum/SED of the NGC 6240 southern 
nucleus outside this range is different from those of M 82 and NGC 253 in three respects 
(Figure 6): (1) the NGC 6240 SED rises more sharply in the mid-infrared; (2) the NGC 6240 
spectrum drops more slowly at < 7 ^m; and (3) the NGC 6240 SED has more luminosity in 
the far-infrared. The third point is less certain because this depends on the far-infrared SED 
of the southern nucleus, which we cannot measure directly. In Figure 6, we estimated the 
southern nucleus far- infrared SED of NGC 6240 from the IRAS 60/100 /^m measurements 
assuming a luminosity ratio of 4 between the two nuclei (i.e., the southern nucleus luminosity 
contribution is 80%) . If the true ratio is larger as suggested by the mid-infrared measurements 
(up to ~8), then the discrepancy becomes larger; On the other hand, if the true ratio is 
as small as the radio measurements indicate (down to ~ 2), the disrepancy, at least in 
comparison with NGC 253, would go away. 

It is interesting that between the two SEDs of the starburst galaxies, the SED of 
NGC 253 is significantly closer to that of NGC 6240. NGC 253 is also known to con- 
tain an AGN, detected in the radio (Ulvcstad & Antonucci 1997) and X-ray (Weaver et al. 
2002), and this might explain the similarity. In fact, as we shall show below, it is possible to 
reproduce the SED of NGC 6240 from that of M 82 by adding spectral components produced 
by a heavily dust-obscured AGN. On the other hand, the AGN in NGC 253 is thought to be 
much weaker than that in NGC 6240, and therefore it may not have a significant effect on 
the shape of the infrared SED. In this case, the difference between the SEDs of M 82 and 
NGC 253, though quite significant (Figure 6), should simply be understood as a variation of 
properties among starburst galaxies. Although the situation is not yet clear, the simplest in- 
terpretation would be the latter: that is, there is a significant variation of properties among 
starburst galaxies, and therefore that everything we see in the infrared spectrum/SED of 
NGC 6240 can still be explained by starbursts. 



4.2. The Hidden AGN 

Given the fact that the AGN in NGC 6240 is hidden behind a Compton-thick wall of 

material (Vignati et al. 1999), it is not surprising that there is no clear AGN signature seen in 
the mid-infrarcd. A column of Compton-thick material {Nh > 2 x 10^"^ cm^^) could produce 
a visual extinction of Ay > 1000 mag, which is simply too large to be penetrated in the 
mid-infrared. In fact, Krabbe, Boker, & Maiolino (2001) have shown that compared with 
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Seyfert galaxies, NGC 6240 is an order of magnitude underluminous in the mid-infrared with 
respect to its hard X-ray luminosity, which suggests that the AGN light is highly absorbed 
even in the mid-infrared. 

We have estimated the maximum AGN contribution allowed by the data using a simple 
model. The mid-infrared radiation of the AGN was assumed to be heavily extincted so as to 
erase the AGN spectrum in the mid-infrared, and the absorbed mid-infrared luminosity was 
assumed to be re-emitted in the far-infrared as thermal radiation from dust. The emission 
resulting from this model is shown in Figure 7 and consists of the following three components: 
(1) starburst emission, (2) residual AGN emission after absorption, and (3) re-emitted AGN 
emission. The M 82 SED shown in Figure 6 was used to represent the starburst. A spectrum 
of NGG 1068 (Sturm et al. 2000) reddened by Ay = 70 mag represented the residual AGN 
emission. The NGG 1068 spectrum was extended to longer wavelengths by extrapolating 
the power-law continuum at 18-45 //m. This extrapolation does not affect the modeling 
results as long as it is much less luminous than the other components in the far-infrared. 
The re-emitted AGN was represented by a black-body emission with a temperature of 55 K, 
a value that produces a good fit to the observed SED (a similar temperature of 57 K was also 
derived by Klaas et al. (2001)). In this particular model, the luminosity of the starburst in 
the restricted 4.5-120 /xm range is 1.8 x 10^^ Lq, that emitted by the AGN after extinction is 
1.1 X 10^^ Lq and that of the re-emitted AGN 1.5 x 10^-*^ L©, which makes the total emitted 
luminosity ascribable to the AGN to be 2.6 x 10^^ L©. The intrinsic AGN luminosity could 
even be larger if the AGN emission is not isotropic in the infrared. This particular model 
assumes an intrinsic AGN luminosity of 3.3 x 10^^ Lq, which was a free parameter. So, if 
this model prediction is correct, ~20% of the AGN luminosity is not seen in our line of sight 
as either re-emitted or residual AGN energy after absorption. 

In this model, the deficiency of the M 82 SED is compensated by the absorption and 
re-emission of the AGN light. As seen in Figure 7, the AGN-related emission nicely adds 
luminosity to the M 82 SED such that the final SED matches that of the observed SED. 
The relative contribution to the total observed luminosity is 40% from the starbursts and 
60% from the AGN at 4.5-120 ^m. If we take the intrinsic AGN luminosity (3.3 x 10^^ L©) 
instead of the observed one (i.e., residual plus re-emitted), the AGN contribution becomes 
65%. Although it is not difficult to reduce the AGN contribution significantly by adopting 
a different starburst SED (e.g., NGG 253), this particular model indicates that it is also 
possible to bring up the AGN contribution level to 60-65% in NGG 6240 despite the absence 
of any obvious AGN signatures in the mid-infrared spectrum. Given various uncertainties, 
these numbers are compatible with the estimates by Lutz et al. (2003) based on a number of 
AGN-starburst diagnostics including the properties of the mid-infrared fine structure lines 
observed by ISOSWS. The purpose of this modeling, however, is not to promote the view 



- 11 - 



that the AGN dominates the luminosity output of NGC 6240 since the model presented here 
is just one example in which we pushed the AGN luminosity to its maximum despite some 
difficulties. The intention is rather to illustrate how subtle the effects can be when the SED 
of a heavily absorbed AGN (the residual AGN emission after absorption plus re-emitted 
components) is added to that of starbursts. 

Perhaps the most serious flaw of this model is the overproduction of the flux at 13-18 /xm. 
It is difficult to suppress the AGN flux in this wavelength range because it corresponds to 
the local minimum of the extinction curve between the 10 and 18 fira silicate absorption. 
However, given our limited knowledge of the extinction curve in this wavelength range, we do 
not yet consider this problem as fatal. Put in another way, if the extinction curve adopted 
here is correct, this is the wavelength range in which we may be able to detect heavily 
dust-obscured AGNs. 

In fact, the choice of an extinction law affects such modeling substantially. For the 
residual AGN spectrum shown in Figure 7, we used the extinction law by Weingartner & 
Draine (2001) for the Milky Way dust size distribution with total extinction to reddening 
Rv — Av/E{B — V) —3.1 and increased the amount of extinction between 5 and 9 //m to 
fit the extinction law found for the hne of sight to the Galactic center (Lutz et al. 1996). 
This modification was necessary to suppress the AGN light and produce a good fit to the 
observed spectrum. Without this modification, the AGN light would have dominated the 
spectrum at < 7 /im with the assumed level of extinction {Ay = 70 mag) . 

4.3. The Limitations of the AGN-Starbursts Power Diagnostics 

When, as seems to be the case with NGC 6240 (Krabbe, Boker, & Maiolino 2001), the 
AGNs are highly absorbed in the mid-infrared, the mid-infrared diagnostics (e.g., Genzel et 
al. 1998; Laurent et al. 2000; Imanishi & Dudley 2000) lose their sensitivity to AGNs. The 
most reliable AGN-starbursts power diagnostics arc apparently the hard X-ray-bolometric 
luminosity correlation for AGNs and the radio-infrared luminosity correlation for starbursts. 
However, neither of these diagnostic relations can be used to estimate the bolometric lumi- 
nosity of an AGN or starbursts with an accuracy better than a factor of a few. The radio- 
infrared luminosity correlation for starbursts is quite tight over a large luminosity range, 
but even in the most recent and extensive compilation (Yun, Reddy, & Condon 2001), this 
correlation has an intrinsic 1 a scatter of almost a factor of two (0.26 in dex). The scatter is 
larger with the correlation between the AGN hard X-ray and bolometric luminosities. This 
means that when the starburst and AGN luminosities are comparable (i.e., within a factor 
of a few of each other), it is impossible to determine which dominates and by how much. 
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The SED variation between the two archetypical starburst galaxies M 82 and NGC 253 

(Figure 6) underscores this intrinsic uncertainty associated with any AGN-starbursts power 
diagnostics. Despite the good match of the mid-infrared spectra, their SEDs differ by a factor 
of a few in the far-infrared, where most of the luminosity is coming out. This would leave 
room for a substantial contribution from a heavily dust-obscured AGN that emits mostly in 
the far- infrared. 

The AGN-starbursts power diagnostics based on mid-infrared data would likely suffer a 
larger uncertainty compared with the diagnostics using hard X-ray and radio data because of 
the larger effects of extinction involved, which are difficult to correct accurately. As seen in 
Table 2, the flux ratio of the two nuclei in the mid-infrared (7:1) is significantly different from 
the ratios in the hard X-ray (3:1) and radio (2:1). If we take the hard X-ray and radio ratios 
as the luminosity ratios of the AGN and starburst components, respectively, then this means 
that the the bolometric luminosity ratio of the two nuclei must be 2-3 regardless of what the 
relative starburst/ AGN luminosity contributions are. However, for some reason (most likely 
due to a larger extinction in the northern nucleus), this is not the case in the mid-infrared 
out to 24 fim. The far-infrared power diagnostics (e.g., the [C II] 158 fim line luminosity vs. 
far-infrared luminosity) are more promising in principle, and according to these diagnostics, 
NGC 6240 is a starburst galaxy (Luhman et al. 1998, 2003). However, unless this correlation 
is significantly tighter than that of the radio-infrared luminosity correlation, we will again 
be left with the factor of a few uncertainty. 

What may be even more problematic is the assumption that the infrared radiation from 
a ULIRG is isotropic even when it contains a dust-obscured AGN. A natural consequence of 
a thick dusty torus model is that most of the AGN radiation is directed perpendicular to the 
plane of the torus. This means that if we are looking at the torus from an edge-on direction, 
which may be the case for NGC 6240 considering the large extinction derived from the hard 
X-ray observations, most of the AGN luminosity may be invisible to us. In other words, even 
if the observed infrared luminosity is almost all due to starbursts, this would not necessarily 
exclude the existence of an equally (or even more) powerful AGN inferred from the hard 
X-ray observations since the true bolometric luminosity of NGC 6240 may be substantially 
larger than calculated from our line of sight. Observationally, the detection of the flat (i.e., 
unreddened) reflected hard X-ray continuum (Iwasawa & Comastri 1998; Boiler et al. 2003) 
and a high-ionization [O IV] hue (Genzel et al. 1998; Lutz et al. 2003) suggests that the AGN 
is not completely covered up in all directions by Compton-thick material, which is consistent 
with a torus-like geometry. 

For a toroidal model to work, it also requires the existence of some material which 
obscures the AGN in the polar directions of the torus and softens the AGN radiation. 0th- 
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erwise, the various emission lines would show excitation levels higher than observed, which 
are those of a LINER in NGC 6240 (Heckman et al. 1983). There are indications in some 
Seyfert galaxies that AGN radiation is significantly softened before reaching the narrow line 
regions (Alexander et al. 1999, 2000), but it remains to be seen if such a model would work 
in detail to reproduce the observed properties of NGC 6240. 

5. CONCLUSIONS 

Using the Keck Telescopes, we have obtained mid-infrared images and spectra of the 
luminous infrared galaxy NGC 6240 with sub-arcsecond spatial resolutions. The main con- 
clusions are as follows: 

1. The mid-infrared (and therefore presumably the far-infrared) luminosity of NGC 6240 
mostly originates from the two nuclei located at the positions of the southern and 
northern radio/hard X-ray Sources. 

2. The southern nucleus is especially luminous and compact (<200 pc in diameter), emit- 
ting 80-90 % of the mid- infrared luminosity of the two nuclei. 

3. The southern nucleus seems to be powered by starbursts. Its mid-infrared spectrum is 
similar to those of the local starburst galaxies M 82 and NGC 253, and its mid-/far- 
infrared SED is broadly consistent with those of the starburst galaxies. 

4. At the same time, because of the slight mismatch between the NGC 6240 SED and 
those of the starburst galaxies, a significant AGN contribution (up to ~60 %) to the 
bolometric luminosity is also possible if the AGN is heavily obscured and emits mostly 
in the far-infrared. 

5. Due to the intrinsic variation of properties among a given galaxy population, there is 
a scatter in the relationships inherent even with the most reliable and accurate AGN- 
starbursts power diagnostics. When the starburst and AGN luminosities are within 
a factor of a few of each other, this makes it impossible to determine in individual 
galaxies whether an AGN or starbursts dominates the luminosity. 
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Table 1. N6240 mid- infrared continuum flux densities 



Wavelength 


Tint 


TotaP 


South^ 


North'^ 


(/xm) 


(s) 


(mJy) 


(mJy) 


(mJy) 


7.9 


360 


494±12 


275±3 


59±3 


8.8 


360 


243±6 


142±1 


29±1 


10.3 


315 


110±8 


67±2 


7±2 


11.7 


360 


316±4 


201±1 


41±1 


12.5 


360 


480±6 


375±2 


58±2 


17.9 


480 


1098±30 


919±9 


119±9 


24.5 


432 


3263±112 


2734±31 


384±31 



^Measured with a 4" diameter circular beam centered 
between the two nuclei. 

^Measured with a 1" diameter circular beam centered 
on the nucleus. 

Note. — The uncertainties listed in the table arc sta- 
tistical only. Photometric uncertainties in the MIRLIN 
data are ±5% for < 20 and ±10% for > 20 //m. 
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Table 2. N6240 continuum flux ratios 



Band 


Fliix ratio 


Ref 


2-10 keV 


2.8*" 


1 


1.1 fiin 


2.7 


2 


1.6 /Ltm 


3.0 


2 


2.2 /xm 


3.5 


2 


7.9 urn 


4.7 ±0.2 


3 


8.8 fim 


4.9 ±0.2 


3 


10.3 /xm 


> 9.6*^ 


3 


11.7 /um 


4.9 ±0.1 


3 


12.5 /um 


6.5 ±0.2 


3 


17.9 


7.7 ±0.6 


3 


24.5 /xm 


7.1 ±0.6 


3 


1.3 mm 


3.2 ±0.3 


4 


2 cm 


2.5 ±0.5 


5 


3.6 cm 


2.2 ±0.2 


6 


6 cm 


2.3 


7 


20 cm 


2.1 


7 



References. — (1) Komossa et al. 
(2003); (2) Scoville et al. (2000); (3) 
This work; (4) Tacconi et al. (1999); 
(5) Carral, Turner, & Ho (1990); (6) 
Colbert, Wilson, & Bland-Hawthorn 
(1994); (7) Beswick et al. (2001) 

"^Corrected for absorption. 

'^We treat the 3.5 ct detection of 
the northern nucleus at 10.3 jim (Ta- 
ble 1) as an upper limit. 
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Fig. 1. — Mid- infrared images of NGC 6240. North is up and east is left. The observed wave- 
length (in fim) is denoted in each panel. The contour maps were produced from Gaussian- 
smoothed images (Gaussian FWHM = 2 pixels). The lowest contour corresponds to 2a above 
the sky level, which was calculated in the non-smoothed image, and the successive contours 
correspond to a factor of 1.4 increases in surface brightness. The diagonal lines from upper 
left to lower right seen in some images (e.g., at 8.8 //m) are due to the noise pattern of the 
array. 
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Fig. 2. — The 12.5 fim image (top row) compared with the HST/NICMOS images at 2.2 ^m 
(middle row) and 1.1 fim (bottom row) from Scoville et al. (2000). North is up and east 
is left. In each row, the left panel shows the contour map while the right panels show 
the image itself with different gray scales to emphasize the two sources in each nucleus. 
The contour maps were produced from Gaussian-smoothed images (Gaussian FWHM = 2 
pixels) . The lowest contour corresponds to 2a above the sky level, which was calculated in 
the non-smoothed image, and the successive contours correspond to a factor of 1.4 increases 
in surface brightness. The two crosses overlaid on the contour plots indicate the positions 
of the two radio peaks (Gallimore & Beswick 2004) with the assumption that the southern 
nucleus is spatially coincident at radio and mid-infrared wavelengths. The slit position of 
the LWS spectra, which has a position angle of 16° (east of north), is also overlaid. 
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Fig. 3. — The two-dimensional LWS spectrum of NGC 6240 covering 9.9-12.4 /xm in tlie 
restframe: (a) the observed spectrum; (b) the continuum-subtracted spectrum, showing a 
strong 11.3 /im PAH feature in the southern nucleus; (c) the spatial profiles of the continuum 
(dotted hne) and 11.3 //m PAH feature (sohd line) produced by summing the 11.0-11.7 //m 
region. Normalization of both components was achieved by making the integration of the 
continuum peak of the southern nucleus equal to unity. The two horizontal white dashed lines 
in (a) and (b) indicate the positions of the northern (upper) and southern (lower) nucleus, 
respectively although the northern nucleus is too faint to see. 
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Fig. 4. — (a) The LWS spectrum (solid line) and MIRLIN photometric measurements (solid 
circles) of the southern nucleus in NGC 6240; (b) The ISOCAM CVF spectrum by Lutz et al. 
(2003) (dashed line) and MIRLIN photometric measurements with a 4" diameter beam (open 
circles). For comparison, the LWS spectrum of the southern nucleus shown in (a) is over- 
plotted (solid line). The horizontal bars with the MIRLIN points indicate the photometric 
passbands while the vertical bars show a conservative photometric uncertainty of ±10%. 
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Fig. 5. — The mid-infrared spectra of the southern nucleus of NGC 6240 are compared 

with those of the local starburst galaxies M 82 and NGC 253. The thick sohd hne is the 
Keck/LWS mid-infrared spectrum of the southern nucleus while the thick dotted line is the 
ISOCAM CVF spectrum (Lutz et al. 2003) scaled to match the LWS spectrum. Note that 
the ISOCAM spectrum contains both nuclei and the surrounding area. The thin dashed and 
dotted lines show the ISOSWS spectra of M 82 and NGC 253 (Sturm et al. 2000) scaled to 
match the LWS spectrum. 
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Fig. 6. — The mid-/far-infrared SEDs of the southern nucleus of NGC 6240 are compared 
with those of the local starburst galaxies M 82 and NGC 253. The ISOSWS spectra were 
extended to 100 /im using the IRAS 25, 60 and 100 /xm measurements (triangles and squares) 
scaled such that the IRAS and ISOSWS agree at 25 /xm. The 17.9 and 24.5 /xm photometric 
points of the NGC 6240 southern nucleus are from MIRLIN while those at 60 and 100 /im 
are from IRAS with a scaling of xO.8 to account for the southern nucleus contribution to 
the total flux. 
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Fig. 7. — SED model fit for the southern nucleus with both starburst and AGN compo- 
nents. The observed mid-infrared spectra and mid-/far-infrared SED shown in Figure 6 are 

reproduced by three components as discussed in the text: (1) starburst — M 82 spectrum 
(dashed fine; open triangles); (2) the residual AGN emission after absorption — NGC 1068 
spectrum reddened by Ay = 70 mag (dotted line); (3) re-emitted AGN — 55 K black-body 
radiation (dash-dot line). The sum of these three components is shown by the thin solid 
line. The diamonds are the measurements by Kiaas ct al. (2001), and were scaled by xO.8 to 
simulate the fiuxes from the southern nucleus. The open and solid circles are the NGC 6240 
photometric measurements shown in Figure 6. 



